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We have carried out calculations of electron transport through a metal-molecule-metal junction 
with metal nanoclusters taking the part of electrodes. We show that negative differential resistance 
peaks could appear in the current-voltage curves. The peaks arise due to narrow features in the 
electron density of states of the metal clusters. The proposed analysis is based on the ab initio 
computations of the relevant wave functions and energies within the framework of the density 
functional theory using NRLMOL software package. 
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Electron transport through molecular-scale systems 
has been intensively studied in the past two decades [![. 
Largely, the unceasing efforts of the research community 
to further advance these studies are due to important ap- 
plication potentials of single molecules as active elements 
in various nanodevices intended to complement current 
silicon based electronics Q. Among various impor- 
tant properties of the electron transport through metal- 
molecule junctions one may separate out the negative dif- 
ferential resistance (NDR) , that is the decrease of the cur- 
rent / while the bias voltage V across the molecule in- 
creases. The NDR efffect was originally observed in tun- 
neling semiconducting diodes 0, 0|- Later, the NDR was 
viewed in quantum dots[5L|6| and metal-molecule-metal 
junctions (see e.g. Refs. B S fi H EH EI El, El, El)- 

Several possible scenarios are proposed to explain 
the NDR occurence in the electron transport through 
molecules. The NDR could appear due to alignment and 
subsequent misalignment of the Fermi levels of the leads 
with molecular orbitals which happens as the bias voltage 
increases llj, [15[ . This may noticeably modify the cou- 



pling of the molecule to the leads. The variations in the 
coupling strengths could serve as an immediate reason 
for the NDR peaks to appear. Also, such peaks could oc- 
cur as a Coulomb blockade induced effect (BLEU, and/or 
they could originate from conformational changes in the 
molecule [17] and electron interaction with the molecule 
vibrational modes 0, EH • It is likely that different mech- 
anisms could play a major part in the NDR appearance 
in different molecular junctions where it was observed so 
far. However, reviewing the available experimental data, 
one may conclude that the most distinguished NDR fea- 
tures in the current-voltage curves (sharp and narrow 
peaks separated by intervals of extremely low conduc- 
tivity, like those reported in the Ref. [8|) are usually 
attributed to the matching-mismatching of the molecule 
energy levels with those of the leads. Keeping in mind 
that precisely such NDR characteristics are potentially 
valuable for molecular electronics applications, we fur- 
ther concentrate on this mechanism. 



FIG. 1: (Color online) Schematic of the considered system 
including large electron reservoirs labelled S and D, respec- 
tively; metal nanoclusters L and R ; and the molecule placed 
in between. 



Commonly, while studying electron transport through 
molecules, one assumes that electrodes are large enough 
to have a feartureless electron density of states below 
the Fermi energy. In the present work we analyze the 
NDR effect in the electron transport through a molec- 
ular junction where the leads are small metal clusters 
whose electron density of states reveals sharp and dis- 
tinct features. This system provides better opportuni- 
ties to analyze the effects of matching-mismatching of 
the molecule energy levels with those of the "nanoleads" 
in the electron transport characteristics. We show that 
under certain conditions such junctions may show very 
distinguished multiple NDR features in the I — V char- 
acteristics. This demonstrates their potential usefulness 
in nanoelectronics applications. 

To simplify the computational procedure we omit 
from the present consideration effects originating from 
Coulomb interactions of electrons and from molecule vi- 
brations. To maintain a steady supply of electrons tun- 
neling through the junction, we assume that the metal 
clusters ( "nanoleads" ) keep contact with large source and 
drain electron reservoirs as sketched in the Fig. 1. The 
latter, however, are separated from the molecule is such 
a way that electrons cannot directly tunnel from these 
reservoirs to the molecule. In the absence of the ap- 
plied bias voltage the whole system is supposed to be in 
the equilibrium state characterized with the equilibrium 
Fermi energy Ep. We write the effective Hamiltonian 
for the molecule (H e ft) in the usual form [J, EH 2(| 21 1: 
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(1) 



where the term Hm corresponds to the molecule itself, 
and Hl ji describe the coupling of the latter to the metal 
nanoclusters. 

Omitting for a while the molecule coupling to the leads, 
one may introduce the retarded Green's function for the 
molecule. The latter is defined by the matrix equation: 

[(E + i V )S-H M ]G R = i (2) 

where I is the identity matrix, S is the overlap matrix, 

Sij = J i>*{v)^{v)d\ (3) 

rj is an infinitesimal positive parameter (rj — > + ) , E 
is the energy, and V'iiV'j are t ne orbitals included into 
the basis set. 

Then we employ the Dyson equation. This equation 
relates the Green's function of the molecule coupled to 
the leads G R (E) to the Green's function of the single 
molecule G R {E). It reads 0: 



g r (e) = {[gK(e)]- 1 -£(e)Y 



(4) 



Here, the self-energy correction consists of two 

terms describing the effect of two clusters: 



H(E) = E L (E) + E R (E). 



(5) 



For convenience, one may introduce a notation 
[Gq(E)]^ 1 — A(E), which allows to simplify the form 
of the Eq. |4j namely: 



G R (E)=[A(E)-t(E)\ 



(0) 



Now, one must calculate matrix elements Aij(E) and 
Y>ij(E) Assuming that the wavefunctions are orthonor- 
malized molecular orbitals, the matrix A(E) is a diago- 
nal matrix: 



Aij{E) = (E + irj-Ei)S i: 



(7) 



where Ei are the energy eigenstates of the molecule. 

The matrix elements of self-energy corrections have the 
form [l[: 
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here, (3 € L, R, Wn^p are, respectively, the coupling 
strengths between "i" -th molecule state and "fc" -th 
state on the left/right metallic cluster (lead), tk,p are 
energy levels of the corresponding leads, and the param- 
eters (Jk,/3 are the self-energy corrections which originate 
from the coupling of the clusters to the large electron 



reservoirs. Their imaginary parts characterize the width 
of the clusters energy levels. The summation over "fc" 
in the Eq. [8] is carried out over the states of the left /right 
cluster. 

When the bias voltage V is applied across the system 
shown in the Fig. 1, this causes charge redistribution, 
and subsequent changes in the energies Ei and €k,p- In 
consequence, the matrix elements A+j and S^- values 
vary as V changes. This affects the electron transmis- 
sion function T given by the expression: 



T = Tr{t L G R f R G A ) 



(9) 



where f L ^ = — 2ImEi.^, and G A is the advanced 
Green's function of the molecule (G A = {G R )^\ . When 
the voltage V is applied, the electron transmission which 
determines transport properties of the molecular junction 
may significantly depend on its value and polarity. 

The electron tunnel current through the junction could 
be written in the form: 

dET(E, V mol ) [f(E - Ms ) - f(E - hd)] ■ 

(10) 

Here, T(E, V mo i) is the electron transmission, given by 
the Eq. [9l f(E) is the Fermi distribution function for 
the energy E. Chemical potentials fJ-s.D arc attributed 
to the source and drain reservoirs, respectively. They are 
shifted with respect to the equilibrium Fermi energy due 
to the applied voltage: 



Ms 



E F + (l- v)eV; 



fJ-D 



E f 



veV 



(11) 



where e is the electron charge and v is the division 
parameter. 

We remark that the transmission T in the Eq. [10] 
actually depends on the voltage applied across the junc- 
tion (Vmoi) ■ The latter may noticeably differ from the 
external voltage V. The difference originates from both 
charge redistribution inside the molecule and the electro- 
static potential drops between the source/drain reservoirs 
and the metal clusters. The current-voltage characteris- 
tics shapes crucially depend on the electrostatic potential 
profile in the considered system. If the voltage mostly 
drops between the large electron reservoirs and the small 
metal nanoparticles included in the junction (V mo i V) 
one may approximate T(E, V) as T(E, 0) . Then the 
applied bias voltage does not change relative positions 
of the energy levels tk,f3 and Ei, and low-temperature 
characteristics should display step-like shapes. These are 
typical for electron tunneling through molecules (see e.g. 
[l3,[2^|). Current increases as the voltage V increases, 
and the NDR does not appear. On the contrary, when 
no significant voltage drop occures between the electron 
reservoirs and the metal nanoclusters (V mo i ~ V), the 
effect of the bias voltage on the energy levels of the clus- 
ters £k,/3 could considerably differ from its effect on the 
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FIG. 2: (Color online) Left panel: The calculated electron 
transmission function through the molecule T(E, 0) within a 
certain range of energies E. Right panel: The calculated cur- 
rent voltage characteristics for the electron tunneling through 
the junction. The curves are plotted asuming T = 50-ftT, and 
V mo i = 0.1 V (solid line) and V mo i = 0.9E (dashed line). 



molecular energy E{. Due to these differencies in the 
voltage induced shifts, the relative positions of the en- 
ergy levels on the clusters and these on the molecule 
vary as the voltage changes. This creates opportunities 
for alignment /disalignment of the molecule orbitals with 
those associated with the metal clusters. Therefore, in 
such a case one may expect the NDR to occur. 

To confirm the above suppositions we carried out cal- 
culations of the tunnel electric current through a junc- 
tion which consists of two copper nanoclusters and an 
alaphatic-saturated hydrocarbon chain situated in be- 
tween them (see Fig. 2). The relevant eigenenergies 
and matrix elements included in the Eqs. [3 E1 ar e com- 
puted using the NRLMOL software package [24j. First 
we performed full self-consistent calculations on the con- 
sidered nanosystem. In addition to solving the Kohn- 
Sham equations self consistently we have optimized the 



geometry as well. Then we use Lowden's method [25 1 
of symmetrical orthogonormalization to construct atom 
centered wannier-like functions [2r3 | from the nonorthog- 
onal atom-centered gaussian orbitals. We reconstructed 
the Hamiltonian matrix in this basis and block diagonal- 
ized it separating out the blocks coprresponding to the 
hydrocarbon chain and the copper clusters. Small off- 
diagonal elements between the L-M and M-R blocks of 
the Hamiltonian describe the coupling of the molecule 
chain to the nanoleads. 

The final matrix was used to determine the various 
matrices needed in Eq. [TlfTTl In the following calcu- 
lations of the tunnel current we assumed V mo i — 0.1 V 
and V mo i = 0.9V, respectively. The resulting I — V 
curves are presented in the Fig. 2. One can see that 
when the electrostatic potential mostly drops between 
the sourse/drain reservoirs and the nearby copper clus- 
ters {Vmoi = 0.1V), the corresponding I — V curve 
reveals a stepwise profile without NDR features. On 



the contrary, distinct NDR peaks emerge provided that 
there exists a significant drop in the electrostatic poten- 
tial between the copper clusters and the molecule tips 
(Vmoi = 0.9V). Also, one may notice that some NDR 
peaks are rather sharp and narrow, and the current peak 
values are much greater than in the valleys between the 
peaks. This is consistent with the experimental data re- 
ported in the §|. 

In conclusion, we have considered a metal-molecule- 
metal tunneling junction where metal leads are nanopar- 
ticles. Due to the extremely small size of these metal 
clusters their electron density of states reveals distinct 
features which are washed out for larger electrodes. We 
did show that NDR peaks could appear due to alignment 
and subsequent disalignment of the energy levels of the 
metal nanoclusters with those of the molecule when the 
voltage applied across the junction varies. The proposed 
analysis is based on simple assumptions concerning the 
electrostatic potential distribution inside the junction. 
It seems unlikely that these simplifications could cause 
qualitative distortions in the considered NDR manifista- 
tions. However, to quantitatively analyze the effect one 
must properly compute the electrostatic potential pro- 
file employing a self-consistent computational procedure. 
Also, one should take into account the effects of molecu- 
lar vibrations which could result in extra NDR features 
superimposed upon those presently analyzed. Neverthe- 
less, we do believe that we showed that the tunneling 
junctions including a molecule placed in between metal- 
lic nanoparticles could exhibit distinct NDR peaks in the 
I — V curves. This makes such junctions useful is de- 
signing nanoelectronic devices. 
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